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ABSTRACT 

We model the hydrodynamic evolution of the plasma confined in a coronal 
loop, 30000 km long, subject to the heating of nanoflares due to intermittent 
magnetic dissipative events in the MHD turbulence produced by loop footpoint 
motions. We use the time-dependent distribution of energy dissipation along the 
loop obtained from a hybrid shell model, occurring for a magnetic field of about 
10 G in corona; the relevant heating per unit volume along the loop is used in 
the Palermo- Harvard loop plasma hydrodynamic model. We describe the results 
focussing on the effects produced by the most intense heat pulses, which lead to 
loop temperatures between 1 and 1.5 MK. 

Subject headings: Sun: activity — Sun: corona 

1. Introduction 

Nanoflares (Parker 1988) are among the best candidates to explain the heating of the 
solar corona, and, in particular, of the coronal loops (e.g. Peres et al. 1993, Cargill 1993, 
Kopp & Poletto 1993, Shimizu 1995, Judge et al. 1998, Mitra-Kraev & Benz 2001, Katsukawa 
& Tsuneta 2001, Warren et al. 2002, 2003, Spadaro et al. 2003, Cargill & Klimchuk 1997, 
2004, Miiller et al. 2004, Testa et al. 2004). 

Although the evidence of nanoflares appears to be well established, it is still unclear 
whether, and to what extent, they really can provide enough energy to heat the whole 
corona (e.g. Aschwanden 1999). More recently, models of nanoflares with a prescribed 
random time distribution of the pulses deposited at the footpoints of multi-stranded loops 
have been proposed (Warren et al. 2002, Warren et al. 2003), and have been shown to 
describe several observed features. 
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According to some models, nanoflares are the result of dissipation in an MHD tur- 
bulence, generated inside closed magnetic structures in the corona, and due to nonlinear 
interactions among fluctuations generated by photospheric motions. Possible evidence of 
turbulent motions has been detected from line broadenings in coronal loops (Saba & Strong 
1991). Most of these models include direct numerical solution of MHD equations in two or 
three dimensions (Einaudi et al., 1996; Hendrix & Van Hoven, 1996; Dmitruk & Gomez, 
1997; Dmitruk et al, 1998; Dmitruk & Gomez, 1999; Buchlin et al, 2003) using relatively 
low Reynolds/Lundquist numbers. Recently Nigro et al. 2004 (hereafter NMCV04) have 
related coronal nanoflares to intermittent dissipative events in the MHD turbulence pro- 
duced in a coronal magnetic structure by footpoint motions. The injected energy is stored 
in the loop up to significant levels in the form of magnetic and velocity fluctuations and 
released intermittently through nonlinear interactions which process these fluctuations and 
generate cascades toward smaller scales where energy is dissipated. The derived probability 
distribution functions of the peak maximum power, peak duration time, energy dissipated 
in a burst and waiting time between bursts are in good agreement with those obtained from 
the analysis of coronal impulsive events (Datlowe et al. 1974, Lin et al. 1984, Dennis 1985, 
Crosby et al. 1993, Shimizu & Tsuneta 1997, Krucker & Benz 1998, Boffetta et al. 1999, 
Parnell & Jupp 2000, Aschwanden et al. 2000a, b). This heating model does not need any 
ad hoc hypothesis, once the loop length and the characteristic Alfven speed, i.e. the strength 
of the ambient magnetic field (if the density does not change much), are fixed. 

In the present work we model the plasma confined in a coronal loop heated according 
to the events dissipation rate and distribution described in NMCV04. We will compute the 
evolution of the distributions of the density, temperature and velocity of the loop plasma 
by means of the time-dependent thermo-hydrodynamic Palermo- Harvard (Peres et al. 1982, 
Betta et al. 1997) loop model assuming the output of the hybrid shell model illustrated in 
NMCV04 as the basis of the heating function. 

In Section 2 we describe the set up of the loop model with the MHD-turbulence dissi- 
pation rate as input heating; in Sec. 3 we show relevant results and discuss them in Sec. 4. 

2. The loop model 

Our purpose here is to model the evolution of the plasma confined in a coronal loop 
under the effect of the energy dissipation predicted in NMCV04. According to their settings, 
we model a magnetic loop, with a total length of 30,000 km. The plasma is described as a 
compressible fluid moving and transporting energy only along the magnetic field lines, i.e. 
along the loop itself. Thus, the magnetic field has only the role of confining the plasma. The 
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loop model assumes constant loop cross-section. 

We use the Palermo-Harvard code (Peres et al. 1982, Betta et al. 1997), a 1-D hydro- 
dynamic code that consistently solves the time-dependent density, momentum and energy 
equations for the plasma confined by the magnetic field: 

dn dv . . 

dv dp d dv 
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with p and e defined by: 



p=(l + (3)nK B T e =^p + nf3 X , (4) 



where n is the hydrogen number density, s the spatial coordinate along the loop, v the 
plasma velocity, mn the mass of hydrogen atom, \x the effective plasma viscosity, P(T) the 
radiative losses function per unit emission measure, f3 the fractional ionization, i.e. n e /nH, k 
the thermal conductivity Spitzer (1962), Kb the Boltzmann constant, and x the hydrogen 
ionization potential. H(s,t) is a function of both space and time which describes the heat 
input in the loop. This function will be described in detail in Sec. 2.1. The numerical 
code uses an adaptive spatial grid to follow adequately the evolving profiles of the physical 
quantities, which can vary dramatically in the transition region and under the effect of the 
evolution. The loop is not symmetric, the apex is at half the numerical grid and there is a 
chromosphere on each side. The boundary conditions at the loop footpoints are the same as 
in Reale et al. (2000). 



2.1. The heating function 

The original version of the Palermo-Harvard hydrodynamic code includes a space- and 
time-dependent heating function, which describes the input of external energy triggering 
transient events (Peres et al. 1987). Several formulations are possible and the code can be 
easily adapted. For the present work, the heating function is given by the output dissipation 
rate of NMCV04 (in the form of a numerical table). 

The model developed in NMCV04 has been derived within the framework of the Re- 
duced Magnetohydrodynamics (RMHD) (Strauss 1976, Zank & Matthaeus, 1992), with the 
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Fig. 1. — Evolution of the average heating rate per unit volume released in the loop. The 
vertical dashed lines mark the times illustrated in detail in Fig. 2. 
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assumptions that: (i) the plasma is permeated by a strong uniform magnetic field B in the 
longitudinal direction; (ii) there is low thermal to magnetic pressure ratio (3p = 8np/B 2 <C 1; 

(iii) the longitudinal scale /|| of transverse velocity v±_ and magnetic field B± fluctuations is 
much larger than the transverse scale l±; indeed, the MHD turbulence is anisotropic (e.g., 
Carbone & Veltri, 1990), the energy cascade being more efficient perpendicularly to Bo- 

(iv) Small amplitude perturbations B±/B = v±/cao < l±/l\\ 1, where cao is the back- 
ground Alfven velocity, commonly assumed of the order of Cao ~ 10 8 cm s -1 , while the 
fluctuating velocity can be estimated using nonthermal broadening of coronal spectral lines: 
v _i_ ~ 3 x 10 6 — 1.5 x 10 7 cm s -1 . Under the above assumptions the set of the RMHD 
equations can be derived; they describe the evolution of magnetic and velocity fluctuations 
in terms of two distinct effects: (a) wave propagation in the longitudinal direction, at the 
Alfven velocity; (b) nonlinear couplings, which generate a turbulent cascade perpendicularly 
to B . The model proposed by NMCV04 (hybrid shell model) includes both these dynam- 
ical mechanisms, but nonlinear effects are described in a simplified way by using a shell 
technique (Boffetta et al., 1999): a Fourier expansion is carried out in the perpendicular 
directions and the resulting spectral space is divided into concentric shells of exponentially 
increasing radius. In each shell velocity and magnetic field fluctuations are represented by 
complex scalar quantities. Nonlinear effects are reproduced by quadratic terms representing 
the interactions between nearest and next nearest neighbor shells; the coefficients are chosen 
so as to conserve 2D quadratic invariants: total energy, cross helicity and squared magnetic 
potential. The equation of the hybrid shell model is written as: 
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where Z°(s,t) = v n ±(s,t) + crb n j_(s,t) (with n — 0, 1, ...,n max and a = ±1) are the Elsasser 
variables; k n = k 2 n the transverse wavenumber, with k = 2tt(L/L_ l ); x = ^/{caoL), where 
the magnetic diffusivity A has been assumed equal to the transverse kinematic viscosity; the 
asterisk means complex conjugate. Lengths are normalized to the loop length L, and time 
to the Alfven transit time £a = L/cao] the velocity v n j_ and magnetic field b n j_ fluctuations 
are normalized to cao and B , respectively. 

The shell technique allows us to describe the turbulence at high Reynolds/Lundquist 
numbers with a relatively small number of degrees of freedom. In particular, we used a 



- 6- 



number of shells n max = 11, with a very small dissipation coefficient x — 10~ 7 . Since 
the longitudinal spatial dependence is retained, the hybrid shell model can describe effects 
of longitudinal resonance. Moreover, it is possible to implement boundary conditions to 
describe the effects of transverse motions at the loop bases. In particular, the system is 
excited through the boundary at s = 0, by imposing a given velocity perturbation at large 
transverse scales, simulating photospheric motions. This boundary perturbation amounts to 
~ 10 5 cm s _1 , is gaussian distributed and has a correlation time t c = 300 s. At the other 
boundary s = 1 total reflection conditions are imposed. The equations (5) are numerically 
solved using second order finite difference schemes, both in space and in time. 

During the evolution fluctuating energy enters or exits the driven boundary, so the total 
energy content in the loop fluctuates erratically in time. At the same time nonlinear effects 
transfer energy to smaller transverse scales, thus building a turbulence spectrum. Dissipation 
takes place mainly at the smallest scales. Occasionally, the velocity imposed at the lower 
boundary drives the loop near to one longitudinal resonance: then, the velocity fluctuations 
increase at the driven large scale shells, enhancing the energy cascade process towards small 
dissipative scales. This process results in a spike of dissipated energy, converted to heat. 
The dissipated power at time t and position s along the loop is calculated as: 



and is the heating input in the loop plasma model (Eq. 3). The hybrid shell model yields 
the energy distribution along the loop integrated in the transverse direction, and provides 
therefore the heat input for the one-dimensional loop model. The power in the whole loop 
is: 



The profile of W(t) contains a sequence of spikes of different amplitudes and durations. The 
space and time profile of the heating function results from the interplay between the external 
driver (photospheric motions), the loop resonance and the nonlinear turbulent cascade. 

The heat spatial distribution is sampled every 0.1 Alfven time. For an Alfven speed of 
2 x 10 8 cm/s, one Alfven transit time is 15 s (NMCV04). The numerical table yields the heat 
distribution per unit time and volume along the loop (sampled every 37.5 km) and span a 
total time of 307.5 ks, i.e. 3.56 days. We assume a circular cross section and an aspect ratio 
d/L=0.2, where d is the cross-section diameter; the cross-section area is A = 2.83 x 10 17 
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cm 2 . Fig. 1 shows a few selected segments of the evolution of the average loop heating rate 
W(t)/(A L); they are essentially zooms of the dissipation power shown in Fig. 1 in NMCV04. 
The heating per unit volume is negligible in the first 1000 s. After this (relatively short) 
transient, the heating is steadily above 10~ 6 erg cm~ 3 s _1 . The evolution of the average 
heating rate is highly irregular, with sharp pulses whose duration spans all time scales from 
few seconds to a few ks. Some pulses resemble flares. Also the pulses intensity is highly 
irregular. Most of them are entirely below 10~ 4 erg cm -3 s _1 . A few of them are higher 
(although mostly below 10 -3 erg cm -3 s" 1 ); in fact, eleven heating pulses reach values well 
above 3 x 10" 4 erg cm" 3 s~ x and occur around 10.5, 22, 25, 57.5, 69.5, 78, 90, 99, 121, 182, 
249 ks, as shown in Fig. 1. The most intense pulse is the seventh one (90 ks) and is higher 
than 10 -3 erg cm -3 s -1 . The high pulses are noticeably less frequent in the second half of 
the heating time interval: nine of them occur in the first 150 ks. Most of these pulses last 
~ 0.3 — 1 ks and are rather peaked. 

The heating rate per unit volume averaged over the whole heating duration is ~ 3 x 10~ 5 
erg cm -3 s" 1 . According to the loop scaling laws (Rosner et al. 1978), for the prescribed 
length this is the heating rate (per unit volume) of a loop at an equilibrium base pressure of 
~ 0.025 dyne cm -2 and a maximum temperature of ~ 5 x 10 5 K. 

Fig. 2 shows distributions of the heating rate per unit volume along the loop sampled 
during the fourth segment in Fig. 1 (from 22.5 to 27 ks, hereafter segment Refl). For 
each time, a couple of distributions are shown, one at 1.5 s from the other. The heating 
distribution is quite uniform for low heating. During the high intensity phase of the heating, 
the distribution becomes less uniform, with large peaks propagating back and forth along 
the loop and extending over ~ 1/5 of the loop. 

2.2. The initial conditions 

Since our scope is to investigate the structure, stability, and observable properties of the 
simulated loop both in time and on the average, the initial conditions ought to be moderately 
important: we should start with an initially cool and empty loop, thereafter entirely governed 
by the new time-dependent heating. For technical reasons, our choice has been to set up 
this condition by letting an initially hotter loop relax to a much cooler condition. The initial 
loop is obtained from the model of Serio et al. (1981) with a uniform steady heating and a 
base pressure 0.03 dyne cm -2 , corresponding to a loop maximum temperature of ~ 5 x 10 5 
K, i.e. the expected average condition of the nano-flare heated loop. In order to let this 
loop relax, we made a preliminary time-dependent simulation assuming zero coronal heating 
in the loop (but keeping the chromospheric heating on, to have stable footpoints). The 
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Fig. 2. — Spatial distributions of the heating rate per unit volume along the loop sampled 
during the fourth segment of Fig. 1, at the times marked in the fourth panel of Fig. 1. The 
dashed lines are the distributions after 1.5 s from the closest solid line. 
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simulation followed the loop evolution for 2000 s, i.e. approximately 2.5 loop thermal decay 
times (Serio et al. 1991). At the end of the simulation, the loop maximum temperature 
decreased to ~ 60,000 K, and the pressure to ~ 1.5 x 10~ 4 dyne cm -2 . A residual velocity 
field was present in the loop, with speeds not larger than 6 km/s, an amply subsonic (Mach 
0.2) value. We took this final status as the initial condition for the simulations with the 
nanoflare heating. 

3. Results 

Our main purpose here is to explore how the dissipation rate described in NMCV04 
can bring a loop to coronal conditions and maintain it. In this perspective we will describe 
in detail the solution obtained in a segment containing a heat pulse of medium intensity, 
specifically the fourth segment (named Refl, between 22.5 and 26.3 ks) in Fig. 1. We will 
also discuss the segment including the highest heat pulse, i.e. the eighth segment (which we 
will call RefH). The solutions in the other segments do not differ much from those that we 
are going to illustrate. 

3.1. Medium pulse 

Fig. 3 shows the evolution of the temperature, particle density, pressure and velocity 
distributions along the loop obtained from the loop simulations during segment Refl. The 
temperature is steadily below 0.2 MK until the pulse at f « 24.5 ks. Then it gradually 
increases due to the enhanced heating. Fig. 3 clearly shows that the effects of the spatial 
heating structure (Fig. 2) are smoothed by the efficient thermal conduction. The pulse 
drives also plasma evaporation from the chromosphere, visible in the density, pressure and 
velocity distributions (the negative velocity peaks indicate plasma moving upwards from the 
far footpoint). The density distributions shows more significant fluctuations traveling along 
the loop. 

For more quantitative information, Fig. 4 shows selected distributions of temperature, 
particle density, velocity, pressure along the loop around the times marked in Fig. 2. Each 
column of the figure shows the distributions along the loop at the exact time, as well as 
100 s before and after this time. In the low heating state (left column), the temperature 
is steadily between 0.2 and 0.3 MK along most of the loop with a profile very similar to 
that of a static loop. Also the density does not change much along the loop and is always 
below 10 8 cm -3 in most of the loop. The distribution of plasma velocity shows fluctuations 




Fig. 3. — Evolution of the distributions of temperature, particle density, velocity, pressure 
along the loop during segment Refl (the fourth one in Fig. 1). 
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with amplitude ~ 10 km/s propagating back and forth along the loop. During the heat 
pulse, the temperature increases to about 1 MK (in ~ 100 s). The distribution at the time 
of the temperature maximum appears to be more peaked than in the cool state and the 
position of the maximum slightly oscillates around the loop apex. At later times (t > 25 
ks), the temperature slowly decreases and its distribution flattens (right panel). Asymmetric 
fronts of plasma evaporation develop as the heating increases (center panel, solid line) and 
the density starts to increase. The density continues to increase even after the temperature 
maximum (right panel), staying above 2 x 10 8 cm" 3 for a long time. During the heat pulse, 
the plasma evaporation fronts are clearly visible also in the velocity profiles: two similar 
strong fronts rise from both footpoints after t=24.8 ks, reaching a speed of about 50 km/s 
at intermediate positions along the loop. Then the plasma noticeably becomes less dynamic. 
During the heating decay, the loop slowly returns to a cool average state around 0.4 MK. 
The plasma velocity continues to decrease until the plasma becomes practically static around 
t=25.5 ks. Then the velocity distribution gets inverted: plasma begins to drain along the 
loop, at very low speed (lower than 10 km/s). The pressure distribution along the loop is 
quite stable in the cool state. When the heating increases, the pressure increases as well 
(together with the temperature and the density). The pressure distribution then settles to 
a very flat distribution during the pulse decay at about 0.04 dyne cm" 2 . 

Fig. 5 shows the evolution of the loop maximum temperature, the loop minimum density 
and pressure, and of the maximum velocity. The first three quantities are typical of the upper 
region of the loop, close to the apex, the last midway between the apex and the footpoint 
of the loop. The evolution of the loop maximum temperature is globally similar to that of 
the average heating (Fig. 1), but much less noisy. Consequently, it is similar also to the 
evolution of the maximum temperature expected from the evolution of the average loop 
heating through the loop scaling laws (Rosner et al. 1978). The former temperature is 
slightly higher (~ 10 %) and decays more slowly than the latter one. The peak temperature 
is different because scaling laws assume a constant and uniform heating, while the actual 
heating function in the simulation is variable and non-uniform along the loop. The slower 
decay is due to the fact that the plasma response to heating decrease is not instantaneous, 
and the cooling processes have their own characteristic times. The density enhancement 
due to the heat pulse of this segment is significantly delayed (~ 300 s) with respect to the 
temperature increase, as typical of loop plasma evaporation. For comparison, Fig. 5 shows 
the equilibrium loop density values as expected from the loop scaling laws. The comparison 
clearly shows the delay mentioned above, but emphasizes as well that during the pulse rise 
the loop is significantly underdense, and becomes overdense in the later decay phase. This 
is expected in dynamically heated loops: while the heating is on, the loop is filling with 
plasma and therefore below the density equilibrium conditions; when the heating stops, the 
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Fig. 4. — Distributions of temperature, particle density, velocity, pressure along the loop 
sampled during segment Refl (the fourth one in Fig. 1) at the three times marked in Fig. 2 
(one for each column). We show the distributions at the time (solid lines), and 100 s before 
(dotted lines) and after (dashed lines). 
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loop cools down but the plasma drains even more slowly. The maximum pressure has an 
evolution in between that of the density and of the temperature, and explains why the plasma 
dynamics is time-shifted with respect to the plasma thermal evolution. Fig. 5 shows that 
the plasma velocity is constantly below 20 km/s except during the heat pulse, when it grows 
to about 50 km/s. These values are well subsonic. 

From the output results of the hydrodynamic simulations, i.e. distributions of temper- 
ature, density and velocity along the loop sampled at regular time intervals, it is possible to 
compute the UV and X-ray emission from the confined plasma. Fig. 6 shows the emission 
along the loop in three representative XUV lines, i.e. Ca X 558 A, Mg IX 368 A, Mg X 
625 A, peaking at logT = 5.9,6.0 and 6.1, respectively, at the same times as the distribu- 
tions shown in the left two columns in Fig. 4. Since the line emission is sensitive both to 
the temperature and to the square of the density, the emission distributions are less uniform 
and fluctuate more. This may be a distinctive signature of this model in loop observations. 
In these lines the loop is visible for a limited time during this segment. In the hottest line 
(Mg X 625 A) it decays very rapidly. 



3.2. High pulse 

In the course of the whole sequence of heating evolution, the most intense heat pulse 
- which we will label RefH - occurs little after time t=90 ks (eighth panel in Fig. 1). 
Fig. 7 shows the evolution of the loop maximum temperature, the loop minimum density 
and pressure, and of the maximum velocity, to be compared with the evolution obtained in 
segment Refl (Fig. 5). The loop maximum temperature reaches 1.5 MK around time t=90.5 
ks. Then it decays below 1 MK, but stays above 0.5 MK for the rest of the segment because 
of the occurrence of other minor heat pulses. The density at the apex reaches about 4 x 10 8 
cm -3 and a pressure of 0.1 dyne cm -2 around time t=91 ks, about 500 s later than the 
temperature peak. The velocity gets above 60 km/s, always amply subsonic. 

Fig. 8 shows the light curves integrated along the whole loop during segment RefH in the 
171 A and in the 195 filter bands of the Transition Region and Coronal Explorer (TRACE, 
Handy et al. 1999). The light curve in the 171 A filter band resembles the evolution of 
the heat pulses (although much smoother). In the 195 A filter band, only the first pulse is 
significant, and only in its initial phase the emission is significant, giving the impression of an 
anticipated evolution. This evolution resembles more closely the evolution of the maximum 
temperature shown in Fig. 7. 
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Fig. 5. — Evolution of the loop maximum temperature, minimum density, minimum pressure 
and maximum velocity along the loop during segment Refl. The dashed lines indicate the loci 
of the equilibrium conditions of the loop according to the loop scaling laws and corresponding 
to the heating evolution in the Fig. 1. In the velocity plot, the dotted line is the sound speed 
at the loop maximum temperature. 
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Fig. 6. — Emission distributions (erg cm 3 s 1 ) along the loop in three relevant XUV lines 
(Ca X 558 A, Mg IX 368 A, Mg X 625 A) during segment Refl at the same times as the left 
two columns in Fig. 4. For the chosen loop parameters 1CT 10 erg cm -3 s" 1 is a reasonable 
threshold for detection. 
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Fig. 7. — Evolution of the loop maximum temperature, minimum density, minimum pressure 
and maximum velocity along the loop during segment RefH. The dashed and dotted lines as 
in Fig. 5. 
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Fig. 8. — Light curves integrated along the whole loop during segment RefH in the 171 A 
(solid line) and in the 195 (dashed line) filter bands of the TRACE telescope. The 195 A 
emission is multiplied by 4. 
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4. Discussion and conclusions 

This work is devoted to exploring the effect of nanoflares due to the magnetic energy 
dissipation through MHD -turbulence on the dynamic and thermal evolution of the plasma 
in a coronal loop. The parameters considered in NMCV04, i.e. and Alfven speed of 2000 
km/s corresponding to a magnetic field of about 10 G in corona, lead to a loop with a typical 
maximum temperature of 5 x 10 5 K. Since coronal loops are typically observed at higher 
temperatures, > 1 MK, here we focus on the effects produced by the most intense heat pulses 
predicted in NMCV04. We compute in detail the hydrodynamics and thermodynamics of 
the loop plasma during the pulses and analyze the results. 

Although the spatial distribution of the heating has significant fluctuations traveling 
along the loop and also rapid fluctuations in time, we find that the plasma is not so fast 
to react and smoothes out the fluctuations both in space and in time. We find that, under 
the effect of a medium heat pulse, the loop plasma reaches T ~ 1MK and density ~ 
0.2 x 10 9 cm -3 . The efficient thermal conduction makes the plasma respond promptly to the 
heating deposition but also smooths the heating fluctuations. The plasma rapidly reaches 
the equilibrium temperature (according to the loop scaling laws) and then cools following the 
decay of the heat pulse. The same evolution occurs for a higher heat pulse, which produces 
a higher peak temperature of 1.5 MK and a higher density of 0.5 x 10 9 cm -3 . The density 
(and pressure) of the plasma shows more significant fluctuations traveling along the loop 
but globally responds on longer time scales. The heat pulses do not last long enough to let 
the plasma reach the thermo/hydrostatic equilibrium: the plasma is underdense during the 
heat pulse and overdense after the pulse with respect to thermal equilibrium. This density 
evolution is a consequence of the impulsive heating (Winebarger et al. 2003a, Warren et al. 
2003). The speed of the plasma driven by the heat pulse is relatively small, largely subsonic, 
and speeds of few tens of km/s occur only for very few minutes. The emission distribution in 
relevant spectral lines may be relatively more sensitive to fluctuations due to the turbulent 
heating and may be used to diagnose this model. For the highest heat pulse, our model also 
predicts the light curves in two relevant TRACE filter bands to be "out of phase" one from 
the other. This phase difference is in qualitative agreement with observations (Winebarger 
et al. 2003b) but also predicted by other different loop models (Warren et al. 2003). 

The heating model used here has very few free parameters (essentially the magnetic field 
strength and the loop length) and depends on basic physical effects. The shell model does 
not yield a detailed description of turbulence, and cannot reproduce the energy distribution, 
in the direction transverse to the magnetic field. However, it should be adequate to describe 
the behaviour of the loop integrated in the transverse direction and the detailed energy 
dissipation along the loop, matching the scope of the Palermo-Harvard loop model. 
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A series of questions are opened by this work. First, characterizing features of the 
proposed heating are the disturbances traveling along the loop. We have shown that obser- 
vations in single spectral lines may be sensitive to disturbances in the loop, but detecting 
such effects may not be trivial with present day instruments. Also, one may wonder on 
the effect of changing the magnetic field strength; can a stronger field lead to hotter active 
region loops or even major flares? Even if the heating function may be modified with a 
simple scaling, this question requires anyhow additional detailed loop modeling, since the 
loop plasma evolves non-linearly under the effect of the heating, coupled with the dynamics 
and the cooling processes. 

As further issue to investigate, we note that the heating function is modified by the local 
plasma conditions, e.g. the density stratification and its time variation (the Alfven speed 
depends on the density). Including self-consistently a feedback of the loop plasma conditions 
on the energy dissipation may easily modify some characteristics of the heating function, such 
as the pulse duration, and thus influence the results. Tackling this question requires to couple 
the hybrid MHD turbulence model with the loop time-dependent hydrodynamic model, a 
task planned for future work. 

This first work paves the path to future works along several lines, such as the time 
decomposition analysis of results and the coupling of the heating and loop models, the 
comparison with observations, encompassing the selection (or acquisition) and analysis of 
observations made of long and regularly sampled image sequences. 
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